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Inspiralling black hole-neutron star (BH-NS) binaries emit a complicated gravitational wave sig- 
nature, produced by multiple harmonics sourced by their strong local gravitational held and further 
modulated by the orbital plane's precession. Some features of this complex signal are easily ac- 
cessible to ground-based interferometers (e.g., the rate of change of frequency); others less so (e.g., 
the polarization content); and others unavailable (e.g., features of the signal out of band). For 
this reason, an ambiguity function (a diagnostic of dissimilarity) between two such signals varies 
on many parameter scales and ranges. In this paper, we present a method for computing an ap- 
proximate, effective Fisher matrix from variations in the ambiguity function on physically pertinent 
scales which depend on the relevant signal to noise ratio. As a concrete example, we explore how 
higher harmonics improve parameter measurement accuracy. As previous studies suggest, for our 
hducial BH-NS binaries and for plausible signal amplitudes, we see that higher harmonics at best 
marginally improve our ability to measure parameters. For non-precessing binaries, these Fisher 
matrices separate into intrinsic (mass, spin) and extrinsic (geometrical) parameters; higher harmon- 
ics principally improve our knowledge about the line of sight. For the precessing binaries, the extra 
information provided by higher harmonics is distributed across several parameters. We provide 
concrete estimates for measurement accuracy, using coordinates adapted to the precession cone in 
the detector's sensitive band. 



PACS numbers: 04.30.-w, 04.80.Nn, 95.55.Ym 



I. INTRODUCTION 

Ground based gravitational wave detector networks 
(notably LIGO Q and Virgo @] ) are analyzing results of 
design-sensitivity searches for the signals expected from 
the inspiral and merger of double compact binaries. 0,0] ■ 
For the lowest-mass compact binaries M = mi + < 
16Mq, the response of the detector to a binary merger 
with arbitrary masses, spins, and even eccentricity is well 
understood, particularly given the detectors' limited and 
low sensitive frequency band [5Hl6j . Though this compli- 
cated signal encodes all information about the binary's 
spacetime [l7| , the amount of accessible information de- 
pends on signal strength (or signal-to-noise ratio) [Ls| . 
Strong signals permit high-precision tests of general rel- 
ativity; fainter signals allow high-precision constraints 
on some binary parameters; while very faint, short sig- 
nals may only constrain the binary's mass. Qualitatively 
speaking, we can distinguish two configurations if they 
are separated by contours 1 — P > l/p 2 , where P (de- 
fined in Sec. IIII[) is the (normalized) ambiguity function 
or "overlap" and p is the signal-to-noise ratio (SNR). 

As higher-order corrections and new physics are added 
to our models for gravitational-wave signals, the func- 
tional dependence on various parameters (such as masses, 
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spins and orientation angles) in the model grows in com- 
plexity. On scales ~ l/p 2 , for astrophysically plausible 
p, the ambiguity function is generally smooth. However, 
it may have more complicated fine-scale structure which 
may not be detectable for expected signal strengths. The 
Fisher matrix approach to estimating parameter errors 
is based on differentiating a waveform with respect to its 
parameters. These derivatives are defined in an infinites- 
imal patch of parameter space, and are thus measuring 
fine-scale structure, which could potentially be mislead- 
ing about larger-scale, observable trends. This point is 
illustrated by Fig. [2] where we fit a quadratic through 
the ambiguity function. The ambiguity function changes 
shape, so fits to small (P > 0.999) and large (P > 0.99) 
regions of parameter space would give rather different es- 
timates of posterior widths and thus parameter accuracy. 
Similarly, if the standard Fisher matrix were calculated 
via finite difference, step sizes on these scales would give 
different results, with the smaller step size giving a mis- 
leading estimate of parameter accuracy for a signal of 
expected strength. 

In this paper we propose a simple effective procedure to 
identify relevant scales and parameter correlations, con- 
struct suitable "effective Fisher matrices" , and estimate 
ambiguity functions at low but nontrivial signal to noise 
ratio. To demonstrate this technique, we examine the 
signal from selected black hole-neutron star (BH-NS) bi- 
naries, described in Section [TXJ In this paper we use all 
available knowledge about the (post Newtonian) wave- 
form, adopting a complete model for the adiabatic qua- 
sicircular inspiral of precessing BH-NS binaries. In par- 
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ticular, we employ all available harmonics and amplitude 
corrections, introducing small but non-negligible changes 
to the ambiguity function. In Section IIIII we introduce 
our unconventional effective approach to the local ambi- 
guity function. Using those tools, in Section llVI we con- 
struct and approximate the overlap for signals similar to 
each reference binary. Motivated by parameter estima- 
tion, we provide explicit expressions for the Fisher ma- 
trix, correlations, and marginalized uncertainties for each 
configuration. For our fiducial configurations, higher har- 
monics principally allow us to improve our knowledge of 
the binary orientation, providing fairly little additional 
information about intrinsic parameters for the amplitude 
scales of immediate astrophysical interest. 

Our results are complicated by coordinate-dependent 
effects, notably extreme sensitivity to the reference fre- 
quency at which parameters are specified. We show the 
choice of reference frequency can reduce (or introduce) 
fine-scale structure into the ambiguity function, similar 
to the effect of higher harmonics. Our effective approach 
can partially compensate for ill-chosen coordinates, such 
as the coalescence phase or initial spins. To reduce but 
not completely eliminate these systematic effects, we ex- 
press our results using parameters specified near the peak 
sensitivity of the detector (here, 100 Hz). 



e.g., Vallisneri [2f|. In many cases, including those sin- 
gular values addressed in the text, the singular value cor- 
responds to a bounded parameter (e.g., an angle). The 
singular value simply indicates that parameter cannot 
be measured. In the phenomenological limit described in 
this paper, precisely zero eigenvalues never occur, unless 
a parameter is constrained by symmetry. 

Our goal in this work is to understand the typical shape 
of the posterior p(X\n, Xq) for n a noise realization and 
A, Ao coordinates in the signal space. Using one notion 
of "typical" would produce an average posterior over all 
noise realizations. Such an average posterior, however, 
could be slightly wider than the posterior from any given 
noise realization. Instead, in this work we attempt to 
characterize the typical shape of any one noise realiza- 
tion. To do so, in effect we "transport" each posterior so 
their peak likelihoods lie at the same point in parameter 
space. In practice, our procedure amounts to ignoring 
noise-realization-dcpcndcnt changes to the posterior. 



II. SIMULATIONS 
A. Amplitude-corrected precessing waveform 



A. Context and prior work 

Several studies of gravitational wave detection from 
merging binaries have employed amplitude-corrected 
waveforms and higher harmonics. Investigations of 
space-based interferometers, such as the Laser Interfer- 
ometer Space Antenna (LISA), have historically used 
complete signal models, accounting for both spin and 
precession [19j . As higher harmonics have a small effect, 
however, most previous studies of ground-based interfer- 
ometers have omitted them, emphasizing spin. When 
included, higher harmonics were explored alone for non- 
precessing signals. Higher harmonics can allow detection 
of signals otherwise inaccessible due to the detector's 
limited bandwidth [2(3, [Hi . The relative amplitudes of 
higher harmonics can probe astrophysical mechanisms for 
generating non-circularity [22[. Finally, higher harmon- 
ics (and precession) are well-known to break degeneracies 
and improve sky localization, particularly for LISA [23l — 

Si. 

Several authors have explored the local ambiguity func- 
tion "beyond the Fisher matrix" , including higher-order 
correlation functions [26l.[27j and projection effects due to 
the local shape of the signal manifold (28[. These meth- 
ods still use explicit derivatives of the ambiguity function 
(via explicit derivatives of the signal) to construct their 
series approximations. 

The Fisher matrix is often nearly or exactly singular, 
making inversion numerically challenging. Several au- 
thors have pointed out that a singular value implies an 
unconstrained parameter, limited only by the prior; see, 



In this paper we construct the post Newtonian (pN) 
gravitational wave signal from a BH-NS binary using the 
lalsimulation SpinTaylorT4 code [3lj . which is an im- 
plementation based on the waveforms described in JjJ Q . 
This time-domain code solves for the orbital dynamics of 
an adiabatic, quasicircular inspiralling binary by using 
the so-called TaylorT4 method (see [16[ for an explana- 
tion of this and similar methods) of evolving the orbital 
phase and frequency supplemented with precession equa- 
tions to track the motion of the spins and orbital plane 
[29j . The orbital phase and frequency evolution includes 
non-spinning corrections to 3.5pN order and spin correc- 
tions to 2pN order. The precession equations are given 
to 2pN order. This binary evolution is terminated prior 
to merger, either when it reaches the "minimum energy 
circular orbit" , or when the orbital frequency ceases to 
increase monotonically. 

At each time, the values of the gravitational wave po- 
larizations measured by a distant observer can be con- 
structed from the orbital phase, orbital frequency and 
the orientations of the spins and orbital plane. We 
can construct either the commonly used "restricted" (i.e. 
leading-order) polarizations which contain only the dom- 
inant second harmonic of the orbital phase, or we can 
construct amplitude-corrected polarizations which con- 
tain terms that oscillate at other harmonics of the orbital 
phase (and also higher-order corrections to the second 
harmonic). Expressions for the polarizations valid for 
quasi-circular, precessing binaries are currently known to 
1.5pN order [l5|, [2i| [3(| ■ Throughout this work, when we 
refer to amplitude-corrected waveforms, we mean that we 
use the 1.5pN accurate polarizations. 
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FIG. 1: Coordinate system for the precessing binary. 

The left coordinate corresponds to the conventional GW ra- 
diation frame. (4>nj) is a polar (azimuthal) angle of 
the total angular momentum (J) with respect to the radia- 
tion vector (N). In the left coordinate fijL (oijl) is a polar 
(azimuthal) angle of the orbital angular momentum (L) with 
respect to the total angular momentum (J). In the right co- 
ordinate, N, J, and X2 are coplanar and the shaded region 
indicates the orbital plane. 



B. Simulation coordinates 

For the precessing binaries, LIGO-scale studies have 
been complicated by poor choice of coordinates, asso- 
ciated with the start of the waveform. The waveform 
generation code of the standard LIGO algorithm defines 
all the geometrical parameter values at the initial fre- 
quency (40 Hz for the initial LIGO and 10 Hz for the 
advanced LIGO), and evolve the binary system to get 
the full waveforms. Specifically, the orbit is described 
at some point, including the spin, orbital angular mo- 
mentum vector, and the orbital phase. By contrast, the 
detector is more sensitive to higher frequencies. Allow- 
ing for the decreasing signal strength with frequency [32j , 
the detector is most sensitive to the instantaneous binary 
configuration at 100 Hz (initial LIGO) and 40 Hz (ad- 
vanced LIGO, e.g., see Fig. 2 in [H]). Motivated by this 
fact, we choose the reference frequency (/ ro f), at which 
the instantaneous orientations of the spins and orbital 
plane are defined, to be 100 Hz. 

One significant effect introduced by setting the refer- 
ence frequency at 100 Hz is related to the orbital phase. 
The ambiguity function is dependent on how we choose 
the reference frequency. In appendix B, we describe an 
example and illustrate the significant effects in detail; sec 
also Fig. El 

Furthermore, following Brown et al.[32j|. we define the 
geometrical parameters to be angles between the radia- 
tion vector N, total angular momentum axis J and or- 
bital axis L as in Fig. [TJ For comparison, other conven- 
tions specify some parameters using the line of sight, a 
vector — N pointing from our detector to the binary. 



parameter 


mi 


m,2 


<■ 4> x 


M c 


V 


non-spinning 
aligned-spin 


10 
10 


1.4 
1.4 


tt/4 0.0 0.0 
tt/4 0.0 1.0 


2.994 
2.994 


0.1077 
0.1077 



TABLE I: Fiducial parameters for the non-spinning 
and aligned-spin binaries. We adopt the chirp mass M c 
and symmetric mass ratio rj instead of individual masses. The 
orbital phase is defined at 100 Hz. 

For non-spinning or aligned-spin binaries, the total an- 
gular momentum is parallel to the orbital angular mo- 
mentum and the orbital axis is fixed. In effect, the con- 
ventional radiation frame is equivalent to the geometrical 
frame. 



C. Fiducial simulations and local coordinates 

In the case of a non-spinning binary, the binary is spec- 
ified by 9 parameters. In this work, we choose masses 
(mi, 771%) as intrinsic parameters, the polar (inclination) 
and azimuthal (polarization) angles (i, ip) of the orbital 
axis with respect to the radiation vector and an orbital 
phase <j) as extrinsic parameters. Because we maximize 
the ambiguity function over the polarization, we need 
not take this parameter into explicit account henceforth. 
Remaining parameters are, the distance to the detector, 
sky position (two angles), and the coalescence time. The 
fiducial values of parameters are summarized in Table HI 
Mass components (M = mi + m%) can be expressed by 
the symmetric mass ratio r\ = m\mijM.' 1 and chirp mass 
M c = Mr/ 3 / 5 , we adopt these parameters in this work. 

If the NS spin is assumed to be 0, the aligned-spin 
binary is specified by 10 parameters. 9 parameters are 
the same as the non-spinning case and the additional 
intrinsic parameter is dimensionless BH spin parameter 
X- The fiducial values are also summarized in Table |TJ 

The waveform of the precessing binary can be defined 
by 12 parameters if the NS spin is assumed to be 0. In 
this work, we consider 77, M c , BH spin \, and the opening 
angle /3jl of the precessing cone as intrinsic parameters, 
&JL, &nj, 4>nj, and the orbital phase <p as extrinsic pa- 
rameters. Because we maximize over the polarization an- 
gle ip, the parameter 4>mj is eliminated from further con- 
sideration. Remaining parameters are the distance, sky 
position (two angles), the coalescence time. Throughout 
this paper the units are solar masses (for M c ); radians 
(for angles); or the natural dimensionless units (for 77, \). 

Motivated by [3^], we adopt a challenging reference 
configuration, where the polarization along the line of 
sight oscillates between circularly polarized (L along N) 
and linearly polarized (L perpendicular to TV). Fur- 
thermore, to explore the extent to which higher-order 
harmonics allow measurement of parameters that only 
weakly impact the signal, we consider two possible sets 
of initial conditions for L along its precession cone. The 
fiducial values of the parameters are summarized in Ta- 
ble HU For easel, the orbital axis is perpendicular to the 
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parameter 


mi 


7T12 


X 


Pjl 


Onj 







configuration 


easel 


10 


1.4 


1.0 


tt/4 


tt/4 


0.0 


0.0 


N ±L 


case2 


10 


1.4 


1.0 


tt/4 


tt/4 


7T 


0.0 


N || L 



TABLE II: Fiducial parameters for the precessing bi- 
nary: We adopt the chirp mass M c and symmetric mass ratio 
j) instead of individual masses. All the extrinsic parameters 
are defined at 100 Hz. For the easel, L is perpendicular to 
N, and along for the case2 at 100 Hz. 



radiation vector at 100 Hz, for case2 it is parallel to the 
radiation vector at 100 Hz. All the parameter values are 
the same between both cases except for a.jL- 



D. Fiducial network 

We assume two identical interferometers placed per- 
pendicular to the incident signal, which is the optimal 
sky position of the source. We also assume the two in- 
terferometers are oriented by 7r/4 related to one another, 
giving comparable sensitivity to both polarizations. For 
the incident waveforms, we assume a zero noise limit to 
understand how similar the signals are. While not realis- 
tic, they avoid introducing complexity of the signal due 
to the source sky position. 



III. DISTINGUISHING SIMULATIONS 

A. Ambiguity function 

In this work, we reorganize the two projections of the 
strain tensor h + = e c +h a b/2 and h x = h a f,e Xya i ) /2 into a 
complex function: 



h(t) = h+(t) +ih x (t). 



(1) 



We coherently compare a fiducial signal ho(t, Ao), where 
Ao indicates a fiducial source parameter set, to a nearby 
signal h(t,X), with parameters A, by a complex overlap 



where (ho\h) indicates the conventional overlap of two 
real functions defined by 



(ho\h) 



df 



Sn(f) 



Bc[h (f)h(f)* 



(4) 



In appendix A, we summarize the differences between the 
real and complex overlaps. 

We note that a change of the polarization angle, ip, 
simple causes a rotation of the argument of the complex 
wave strain function, h(ip) — e~ 2% ^h(ip = 0). Thus it is 
trivial to find the value of ip which makes the complex 
overlap purely real, so that 



Im{h \ti) = 0, 

Re(h \h') = (h Qt+ \h' + ) + (h Q>> 



(5) 



\h' 



and the value of the complex overlap maximized over 
polarization angle ip is simply 



max^ Kc(ho\h) = \(ho\h)\. 



(0) 



The complex overlap (like the real-valued overlap) can 
also be maximized over the coalescence time t c via an 
inverse Fourier transform as described in j3^|. In partic- 
ular, one uses the fact that 



h(t c =t) = h(t c = Q)e~ 2mft 



(7) 



and notes that the inverse Fourier transform of the com- 
plex overlap integrand in Eq. @ will compute the com- 
plex overlap for all possible coalescence times of h at once 



(h \h(t c =t)) 



df 



— oc Sn(f) 



[V/)M/)*]e 2 ^ 4 . (8) 



The (normalized) ambiguity function P(Ao, A) between 
two waveforms ho(t, Ao) and h(t, A) is then defined as the 
complex overlap maximized over polarization angle and 
coalescence time, 



P(Aq, A) = max tcj 



\(ho\h)\ 



VPoP) 



(9) 



Unless otherwise noted, all overlaps are maximized in 
time and polarization. This is different from maximizing 
over orbital phase </>; see appendix A and B. 



(h \h) 



df 



S n (f) 



[h (f)HfY 



(2) 



where h(f) is the Fourier transform of h{t) and S n (f) is 
a detector strain noise power spectrum. For simplicity, 
we adopt a semianalytic initial LIGO sensitivity |34|.|35|. 
As pointed out by [33j . this complex- valued expression 
characterizes the ability of a network to distinguish sig- 
nals. The real part of the complex overlap corresponds 
to a linear sum of the conventional real overlaps of the 
two gravitational wave polarizations: 



Re(h \h) = (h ,+ \h+) + (h . x \h x ), 



(3) 



B. Likelihood 

The detector noise N(t) is assumed to be a stationary 
and Gaussian process. Given the detector output S(t) = 
H(t, Ao) + N(t) representing a real- valued signal in real- 
valued noise, the probability for the noise to have some 
realization Ao is I37fl 



(10) 



The posterior probability that the gravitational wave sig- 
nal is characterized by the parameters A, can be ex- 
pressed by p(X\S) oc p(X)L(S\X), where p(X) is the prior 
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probability that the signal is characterized by A, L(S\X) 
is the likelihood, which can be written by 



L{S\X) = C x exp 



-\{S-H{X)\S-H{X)) 



(11) 



where C is a proportional factor which, for simplicity, we 
assume to be 1 in this work. 

Since we consider the complex strain, by choosing the 
appropriate polarization angle we shall write the detector 
output for the detector 1 and 2. 



S 1 =H++N 1 , S 2 = H > 



N 2 



(12) 



also 



s = Si + iS 2 , h = H + + iH x , n = iVi + iN 2 . (13) 

The probability for the noise to have both realizations 
Ni and N 2 is 

p(N l7 N 2 ) oc e -{ifi\Ni)/» e -(N a m/2 (M) 

_ e -Re<JV 1 +iiV 2 |iV 1 +iiV 2 )/2 

Finally, using Eqs. (p2]-[ll]), Eq. (fTT|> can be expressed 
by the complex signals: 



L(s\X) = exp 



1 



Re(s - h(X)\s - h(X)) 



(15) 



Substituting s = ho + tiq into this equation [38j, the 
likelihood is 



exp 
exp 



-Rc(n + ho — h\n + ho - h) 



(16) 



-Re{(/i - h\ha - h) + 2{n Q \h Q - h) 



+(«-oK)} 



where the second and third terms in the a square bracket 
depend on the noise realization. They shift the position 
of the maximum likelihood but only weakly change the 
shape of the likelihood curve. In the limit of high SNR, 
these noise-dependent terms can be neglected, so, 



L 



exp 



exp 



1 



-Re(/io — h\ho — h) 



-^{(ho\ho) + (h\h)-2Re(h \h)} 



(17) 



This equation corresponds to the case where two detec- 
tors are placed to have the maximum response to the 
incident two polarizations [for the detector placement, 
see Section D] . While, Eq. (fTTj) corresponds to one 
detector response to one polarization; see appendix A. 

Using Eqs. @ and ©, and the SNR defined by p 2 = 
(ho\ho) = (h\h), the log likelihood can be expressed by 
our complex overlap convention: 



lni 



V(i-p) 



(18) 



where we assume the same strength for both signals. 

For a given log likelihood, the scale of interest of the 
ambiguity function depends on the signal strength p 2 : 



1 - P < 



(19) 



By approximately identifying the L ~ 1/e surface of the 
likelihood, this condition allows us to estimate the set 
of parameters A which cannot be distinguished from Ao 
with a signal amplitude of p using a signal model and 
noise curve that produces an overlap -P(Aq, A). 1 



C. Fisher matrix 

If A is close to Ao, we can write h - 
in the error AA; = Ao — A 



ha 



dh 

8Xr 



A A, 



h to the first order 



(20) 



So, in the limit of high SNR, the likelihood [Eq. ([X7]l ] is 
given as L = exp(— AA; AAj/2), where is 



. dh 



dh 
dX~ 



(21) 



This definition is analogous to the standard Fisher ma- 
trix, except that it is derived from the complex overlap 
and therefore contains information about both polariza- 
tions. If we assume that the prior p(X) is uniform, the pa- 
rameter estimation errors AAi (i. e., the posterior proba- 
bility density function) can be expressed by the Gaussian 
distribution 



p(AAi) = Ne- T ^ AX ' AX ^ 2 



(22) 



where N = ^/dct(r/27r) is the corresponding normaliza- 
tion factor. 

Using another expression relating the Fisher matrix to 
the log likelihood [26|, |39| and Eq. (fT8"]l , we can write 



F — — 

Ay — 



d 2 lnL(A) 



= P 



dXidXj 
,d 2 (l-P) 



(23) 



dXidXj 



= p 2 T 



i j ■ 



where Ao is the fiducial value of source parameter. Here, 
we define the normalized Fisher matrix Ta. 



1 More properly, the probability F(Lo) = J l>Lq p{x)d\ of having 
likelihood > L lets us create a confidence volume for any target 
confidence level. Because the probability p = F(L ) depends 
sensitively on L, we anticipate the edge of this confidence inter- 
val will depend weakly (e.g., as F~ 1 [p] oc y/\np) on the precise 
probability used to define the threshold. 
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For Gaussian noise and high SNR, the inverse of the 
Fisher matrix is the covariance matrix (E^) of parameter 
errors. The measurement error (c^) of each parameter 
and correlation coefficient (cy) between two parameters 
are defined as 



(24) 



where pi and pij are fitting coefficients and P u 
We calculate the effective Fisher matrix as 



jj 



The correlation coefficients are p-independent but of- 
ten sensitive to small changes in T. Conversely, the mea- 
surement error is inversely proportional to p. For the 
purposes of illustration, we adopt p — 10 whenever we 
calculate er,-. 



D. Relevant scales and effective approach 

The Fisher matrix formally involves derivatives, i.e., 
infinitesimal variations of a parameter dX. In this work, 
we compute an effective Fisher matrix by considering fi- 
nite variations S X on scales which give physically observ- 
able changes to the ambiguity function, P. 

To understand the variability on multiple scales we 
plot a one-dimensional ambiguity function of M c for the 
leading-order amplitude, non-spinning binary in Fig. [2] 
In this figure, the ambiguity function P is calculated via 
Eq. ((9]), changing only M c and fixing all other parame- 
ters to be the same for both signals. For comparison, we 
plot quadratic fits 2 to the ambiguity curve at different 
scales of P > 0.99 and P > 0.999. These are the scales 
of interest for signals with strength p 2 ~ 10 2 or 10 3 (see, 
Eq. (TT9")) '). [Note that Fig. [2] is computed with the refer- 
ence frequency at 40 Hz. For our results, except for this 
figure, we choose the reference frequency at 100 Hz. The 
structure illustrated here is present but less significant at 
100 Hz; see appendix B.] 

The shape of the ambiguity function has structure on 
multiple scales. The neighborhood of P > 0.999 suggests 
a much sharper peak than what is seen at the P > 0.99 
scale. A Fisher matrix computed from formal waveform 
parameter derivatives (defined in the limit dX — > 0) or 
finite difference such that P(Xq, Xq + SX) > 0.999 can be 
overly optimistic about how well A can be measured for 
a signal with p 2 ~ 10 2 . 

Therefore, we wish to define an effective Fisher ma- 
trix from the curvature of the ambiguity function on the 
scales of interest. For example, in the case of two param- 
eters, the fitting function P* is 



( r u)eff 



d 2 P* 
dXidXj 



(26) 



In some cases, especially when a parameter is poorly 
determined, the variation of the ambiguity function with 
a parameter may not be well-described by a quadratic. 
See, for example, the top panel of Fig. [5] Therefore, 
we find it useful to employ an "iterative" procedure to 
find the parameters that are amenable to a quadratic fit. 
For each parameter x, we compute the one-dimensional 
curve P(x, x + 8x) and fit it against 1 — (T xx ) e s8x 2 /2. 
This determines the diagonal elements of the effective 
Fisher matrix, (T xx ) c r . We discard any parameters that 
are poorly fit by this quadratic (checked either "by eye" 
or a with quantitative threshold on the goodness of fit). 
For the well-fit parameters, we determine the off-diagonal 
elements of the Fisher matrix by computing the two- 
dimensional surface P(x, y; x + Sx, y + 5y) and fitting it to 
a function of the form like Eq. (j2"5")) while using the val- 
ues from the one-dimensional fits, i.e. p x = — (T xx ) e fi/2. 
We note this method is used primarily as a sanity check 
to identify any parameters which induce obviously non- 
quadratic variations in P. 

Once we have identified the space of all reasonably 
quadratic parameters, the ambiguity function on that 
space can be approximated as 



P*(A , A + SX) = 1 - (fy-)eff 5X 3 /2 



(27) 



P* = P m ^+Pi5\l+p 2 5X 2 +p 12 5X 1 5X 2 



(25) 



Rather than using the iterative procedure described 
above to find the elements of (Tij) e s one at a time, it 
is straightforward to use a standard least-squares fitting 
technique to simultaneously solve for all of the (Tij) Q ff. 
This will also give a better global approximation to the 
ambiguity function than the iterative approach. As an 
example of the small yet noticeable differences between 
these two procedures, Table IVT1 compares the results for 
computing the effective Fisher matrix from both "itera- 
tive" and "simultaneous" fits to the ambiguity function. 
Everywhere else in this work (Tables HH [GO EO ED EED 
and HH the effect ive Fisher matrix is computed by si- 
multaneously fitting all coefficients. 

In the cases where P is not well-described by a 
quadratic (e.g., see Fig. [5]), we can adopt more compli- 
cated expressions to characterize the functional depen- 
dence of P when these parameters are varied, both in iso- 
lation and in correlation with well-constrained variables. 
As a concrete example, in the absence of higher harmon- 
ics the line of sight from the binary is both weakly con- 
strained by observations and nearly separable in P(Xo, X) 
from other degrees of freedom 3 . Specifically, ignoring 



2 Since the posterior function is a normal distribution in the limit 
of high SNR (see Eq. Q221 ). a quadratic fitting function usually 
best fits the log likelihood function with a fiat prior. 



3 Approximate separability of the line of sight from other degrees 
of freedom follows only in our well-chosen coordinates, where the 
binary configuration is specified at 100Hz. 
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maximization in time and phase, the overlap P(N, N') 
between any two lines of sight can be well-approximated 
by Eq. (Blc) from [33|: 



P. 



angles 



|e 2 *y a (0)Y 2 (0 / ) + e- 2 ^r_ 2 (0)y_ 2 (0')l 



y/(\Y 2 (0)\ 2 + \Y-2(O)\ 2 )(\Y 2 (e>)\< 



\Y-2(e>W) 

(28) 



( 21 

where we use the shorthand Y m = Y 2m (n) and simi- 
larly for to reduce superfluous subscripts and where 
we factor out the common e lm ^ from Y lm 2 ^ . This func- 
tion has wide, nearly flat cxtrcma in n for each fixed n' . 
On the other hand, in the absence of higher harmonics 
the line of sight has little impact on the waveform phase 
versus time away from the orbital plane. We can there- 
fore approximate the ambiguity function for P > 0.99 in 
the top panel of Fig. [5] by 



P — -Panglcs[l - -^{rab)csSX a SXb} ~ T aN SX a SX N , (29) 



where the N index varies over the line-of-sight parame- 
ters (9, </>) and the a, b indices vary over the other param- 
eters and r a jv — for a ^ N. With higher harmonics, 
the functional form above [Eq. ([29]) ] is weakly perturbed 
by additional angular terms of the form 



P(Xq, X) ~ P a nglcs[l ~ ^C^aft) off SX a 5X b ] - T aN 5X a 5X N 



(1 - cos l) G(j)(j) (0 - cj>of 
G cc (cos t — cos tg 



(30) 



- G<t> c {4> - 0o)(cosi - COSi ) 

where G a b is a matrix with G^ ;C ~ 0. This approxima- 
tion both factors out the leading-order angular depen- 
dence and adds additional angular terms with parameter- 
dependent coefficients, designed to correctly reproduce a 
<fi- independent result when cost ~ 0. Although these 
terms allow us to correctly reproduce the non-ellipsoidal 
contours seen in the bottom panel of Fig. [3 Tables NT\ and 
IVIII show that this complicated structure only marginally 
improves the overall fit compared to a purely quadratic 
approximation P ~ 1 — (T a b)esSX a 6Xb/2. Fit parameters 
for this more complicated functional dependence are not 
presented here. 

We also considered a more general fit, treating P max as 
a parameter. While this parameterization has a signifi- 
cant aesthetic advantage - its effective Fisher matrix is 
roughly scale-independent when / rc f ~ 100Hz and agrees 
with analytic calculations - it systematically underesti- 
mates P in the neighborhood of the maximum. As both 
fits work well globally, we favor the simpler procedure 
and adopt P max = 1 except for Table IVl 



E. Comparing to standard Fisher matrix results 

Despite subtle differences associated with time domain 
versus frequency domain waveforms, the complex over- 
lap, higher harmonics, and the line of sight, our results 
for the effective Fisher matrix are directly comparable to 
earlier results calculated with the stationary phase ap- 
proximation pol ]. For example, for emission along the z 
axis, both the real and complex strain have the form 4 
h = i4e" 21 * for / > 0. As a result, the Fisher matrix 
can be well-approximated by an identical average over 
frequency: 



SZo d f\~ h \\d a *)(d b *)/S h 



(31) 



where we neglect derivatives d a A as small compared to 
the leading-order phase dependence. Thus, each compo- 
nent of our Fisher matrix must resemble previous results. 
In fact, as the general definition [Eq. ([2"T])] suggests, each 
component of T a (, (unlike £ = T _1 ) depends only on 
the local response to changing two parameters A a , no 
matter how many parameters exist in the model. There- 
fore, the Fisher matrix for an identical model with more 
parameters will have, as a submatrix, the Fisher matrix 
for the smaller model. By contrast, other methods for ex- 
pressing uncertainty like the covariance £ depend simul- 
taneously on all terms in T. For low-mass binaries, the 
Fisher matrix T is well-known to be poorly conditioned, 
with eigenvalues spanning several orders of magnitude. 
We therefore preferentially compare the component-by- 
component Fisher matrix, rather than the covariances £, 
when comparing results. When presenting results, we 
provide several significant figures to insure all eigenval- 
ues of r remain positive-definite. 

Our complex overlap maximizes over time and polar- 
ization. The analytic Fisher matrix calculated from the 
stationary phase approximation [Eq. (|31[) ] has time and 
phase as parameters. To account for maximizing over 
those parameters, we transform the full Fisher matrix T 
to a smaller-dimensional matrix which projects out those 
dimensions: 



(r a b) r 



1CD 



\Pcd] 



FacQcD^Cb 

-1 



(32) 
(33) 



where C, D run over the t, <f> variables and a, b all other 
variables. In these expressions, the matrix Q is the in- 
verse of the projection of T a b into the t, <fi subspace. 



4 The two differ for / < 0: the complex strain has h(f < 0, z) = 0, 
while the real strain has h(—f) = h(f)*. 
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FIG. 2: One dimensional ambiguity function for the 
non-spinning binary, showing variability on multiple 

scales. Dotted line shows P calculated from Eq. [5] for the 
leading-order waveforms, where we compare our fiducial sim- 
ulation with slightly offset analogs, changing just the chirp 
mass. For comparison, the two solid curves show quadratic 
fits to the ambiguity function, on scales of P > 0.99 (blue) 
and P < 0.999 (red). Note that this result is computed with 
the reference frequency at 40 Hz, a default choice that accen- 
tuates this scale dependence; see appendix B. In this paper we 
adopt a reference frequency at 100 Hz to mitigate the change 
in scale shown here; see Fig. [3] for the comparable result in 
that case. 



F. Comparing to posteriors 

Standard parameter estimation techniques like 
Markov-Chain Monte Carlo produce samples of the 
full posterior probability distributions, including post- 
processed data products like one-dimensional standard 
deviations o~ a and two-dimensional covariances cy. 
For strong signals with well-isolated probability dis- 
tributions, our one-dimensional standard deviations 
and covariances are directly comparable, for identical 
binaries. For fainter signals with broad probability dis- 
tributions, our results will describe part of the posterior, 
in the neighborhood of one extremum. 

For brevity, we have explicitly eliminated two param- 
eters - event time and polarization - and make no pre- 
dictions about any correlation including them. We will 
revisit these parameters, along with asymmetric detector 
response, in a subsequent publication. 



G. Numerical and systematic effects 



FIG. 3: Two methods for fitting the ambiguity func- 
tion: Illustration of differences between our default fitting 
technique (blue curve), which fixes P max = 1, and an alter- 
nate fitting technique that lets P ma x 7^ 1, shown here for the 
P > 0.999 scale (dotted curve). These two methods produce 
comparable global fit functions, but the fitting parameters for 
r differ by tens of percent; see, e.g., Table FVl While numer- 
ical error can produce fluctuations in the overlap (e.g., due 
to insufficient sampling rate; we use 65536 Hz), the change in 
shape shown here is resolved. 

cial small-scale structure; to avoid this effect we sample 
at a variety of data rates, typically either 8192 Hz (for 
P > 0.99) or 16384Hz (for P > 0.999). Finally, the am- 
biguity function can also be impacted by our choice for 
the starting and ending frequency. For our calculations, 
we start integrating the waveform and integrate over all 
power above /start = 40 Hz. In our experience, this pro- 
cedure best mimics the real data processing used in initial 
LIGO searches. However, a not-insignificant amount of 
power is present between 30 and 40 Hz; if included in 
the integral, the overlap would differ by AP ~ 10~ 3 , 
comparable to some fine-scale structures of interest. At 
the other extreme, we terminate our evolution at the 
minimum-energy circular orbit (MECO), where the bi- 
nary energy ceases to decrease monotonically. 

One small, subtle, but important effect is the nonzero 
overlap of the waveforms along the ±z axis 5 . For ex- 
ample, for our non-precessing binary, we find that for 
otherwise identical parameters 

max tc> | (h(z)\h(-z\t c , VO) I si 1.7 x 10" 3 (34) 

Equivalently, in the language of single-detector real over- 
laps, the sine and cosine chirps are not precisely orthog- 
onal, for the same orbital phase 6 . To a first approxima- 



At the very smallest scales, delicate implementation- 
dependent choices can also introduce artificial structure 
into the ambiguity function. We have already extensively 
described how the choice of reference frequency intro- 
duces (coordinate-dependent) structure. Less physically, 
the sampling rate for the waveforms can produce artifi- 



5 Such waveforms would have a real overlap of unity, but the com- 
plex overlap is expected to be zero. See Appendix [XI 

6 The real overlap is usually performed in conjunction with an 
explicit maximization over time and orbital phase, for a single 
harmonic. For that situation, this subtlety does not occur. 
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FIG. 4: Comparison of the ambiguity contours be- 
tween the leading-order (solid line) and higher-order 
(dotted line) waveforms for the non-spinning bina- 
ries. 5X is defined by a difference from the fiducial value in 
the Table [TTJ The lines correspond to P = 0.99 (top panel) 
and P = 0.999 (bottom panel). Higher-order waveforms only 
marginally reduce the ambiguity contours. 



tion the h(±z) signals are basis waveforms; the wave- 
form along any line of sight is a superposition of the two. 
Because these two signals are not orthogonal, the am- 
biguity function generally has fine-scale structure with 
AP ~ 10~ 3 , associated with the overlap of these two 
directions. For example, on this scale and below, the 
overlap between two non-precessing waveforms with just 
I = \m\ = 2 emission is no longer well-described by 
Eq. (|28p . Instead, the ambiguity function gains addi- 
tional fine-scale structure in angle. While extremely use- 



ful, for our purposes this result means that on sufficiently 
small scales 1 — P < 10~ 3 , the complex overlap will have 
additional structure compared to "conventional" inves- 
tigations of single-detector, optimally-oriented overlaps 
(i.e., overlaps of two real h+ signals, extracted along z). 
In particular, this nonzero overlap is partially responsible 
for the small-scale structure seen in Fig. [3J 

Because of the many subtle interpretation and imple- 
mentation issues associated with the smallest ambiguity 
scales, while we investigate the value of effective fitting 
to fine scales (e.g., P > 0.999), for simplicity we empha- 
size results for the scale relevant to most detection events 
(P > 0.99). 

IV. RESULTS 

Using a small set of fiducial simulations, we compare 
non-precessing and precessing signals against their im- 
mediate neighbors, mapping out an ambiguity function 
in each n-dimensional parameter space. 

Higher harmonics perturb the ambiguity function by a 
quantifiable amount (i.e., ST a {, oc v p for p depending on 
a, b and the harmonic) . As has previously been shown 
elsewhere, we find that higher harmonics break degen- 
eracies present in non-precessing, leading-order signals 

[II 11 El . 

As described below, we generally find small but signifi- 
cant scale- dependent disagreement with the conventional 
stationary-phase Fisher matrix calculation, even in the 
absence of higher harmonics. Motivated by Figs. [5]and[31 
as well as the Appendix and Figs. [KflandlTTl we suspect 
that most scale dependence is introduced by sub-optimal 
coordinates and can be minimized by a better choice of 
reference frequency. Despite our best attempts to find 
coordinates well-adapted to the problem, the change in 
r going from P ~ 0.99 to P ~ 0.999 for leading-order 
waveforms is comparable to the change in T going from 
leading-order waveforms to higher harmonics. 

A. Zero spin 

For a system without spin, higher harmonics princi- 
pally provide information about the line of sight. For 
clarity, we will first discuss the most immediately relevant 
scale (P > 0.99). Comparing the solid (without higher 
harmonics) and dotted (with higher harmonics) curves on 
the top panel of Fig. @J we immediately see that higher 
harmonics provide little new information about intrinsic 
parameters, all other things being equal. Equivalently, 
looking at Table IIII1 the effective Fisher matrix V on 
the two-dimensional parameters M c ,rj without and with 
higher harmonics are similar to each other, as well as to a 
standard Fisher matrix computed using stationary phase 
approximation waveforms (labeled as P ~ 1). By con- 
trast, as illustrated by the dramatic difference between 
the top and bottom panel in Fig. [5J higher harmonics 
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JJ1N (JILICI 


leading-order 


higher-order 


fitting scale 


P > 0.99 


T~t A AAA 

-P > 0.999 


P ~ 1 


P > 0.99 


7~> --. A AAA 

P > 0.999 


parameter 


71 ,T ~ 

M c rj 


M c rj 




M c 77 


M c 77 


(r<i)eff 


M c 3012±7 -5505±13 
77 - 10675±25 


3621±69 -6474±125 
12478±228 


2547 -5314 
11954 


3243±8 -6070±13 
12151±27 


3932±77 -7224±143 
14367±270 


Cij 


Af c 1.00 0.971 
n - 1.00 


1.00 0.963 
1.00 


1.00 0.963 
1.00 


1.00 0.967 
1.00 


1.00 0.961 
1.00 




0.00760 0.00404 


0.00618 0.00333 


0.00735 0.00339 


0.00689 0.00356 


0.00578 0.00302 



TABLE III: Effective fitting parameters for a non-spinning binary, at different scales. The fiducial values of 
parameters are (M c , 77) = (2.994, 0.1077). Fitting parameters and uncertainties are calculated from a least-squares fit, treating 
the line of sight as fixed. For the one-dimensional errors ffi, we adopt p = 10. Result for P ~ 1 is calculated by Eq. (|21[) using 
h+ polarization of the SPA waveforms in M c , 77, t c and re f , followed by analytic maximization over t c and 4> Te f ; see the text for 
details. [The effective Fisher matrix corresponds to derivatives of an overlap maximized over t c and ip.] Systematic differences 
exist between the SPA waveforms and real overlap used in our numerical calculations and the analytic result shown for P ~ 1. 
Nonetheless, all methods largely agree: higher harmonics provide fairly little additional information about the chirp mass and 
mass ratio correlations,with all other parameters fixed; see Fig. [3] For comparison, this table provides fitting parameters for P 
for two volumes (P > 0.99 and P > 0.999). Because of the effects described in Fig. [2] Section [ill Gl and Appendix [Bl the two 
fits disagree. 



fitting scale 


P > 0.99 


parameter 




M c 


V 


l 






M c 3883±38 


-7236±66 


-0.08350±0.48 


0.08169±0.23 


(Fij ) c ff 


'/ 




14209 ±138 


0.1181±0.93 


1.172±0.37 




i 






0.02723±0.0009 


-0.001413 ±0.00052 




4> 








0.03196 ±0.0003 




M c 


1.00 


0.976 


-0.00313 


-0.271 


Cij 


'/ 




1.00 


-0.00494 


-0.276 




i 






1.00 


0.0474 












1.00 






0.00740 


0.00387 


0.607 


0.583 



TABLE IV: Effective fitting parameters for a non-spinning binary, at different scales. Higher order: Fitting 
parameters in this table are calculated by least-squares to a quadratic form P ~ 1 — (Iij) e ff <5A l <5A J /2. Measurements of intrinsic 
(M c ,rj) and extrinsic (line of sight) parameters separate: with only one exception (F^ctf, all off-diagonal terms coupling the 
line of sight and intrinsic parameters are consistent with 0. We anticipate a slightly different choice of reference frequency 
will eliminate the small residual correlation that remains. Fitting a general form described in Eq. (|30[) that accounts for the 
manifestly nonquadratic behavior shown in Fig. [S]leads to comparable results: a separable fit (i.e., G c( p = and F a jv — 0) that 
performs little better than the quadratic form used above. 



pN order 


leading-order 


higher-order 


fitting scale 


P > 0.99 


P > 0.999 


P~ 1 


P > 0.99 


P > 0.999 


p 

i max 


0.999767±1.09 x 10~ 5 


0.999884±0.78 x 10~ 5 


n/a 


0.999775±1.01 x 10~ 5 


0.999888±0.86 x 10~ 5 


parameter 


M c 77 


M c 77 


M c 77 


M c 77 


M c 77 


( r «) 8 ff 


M c 2899±8 -5295±14 
77 - 10261±28 


2980±62 -5325±111 
10276±207 


2546 -5313 
11954 


3125±8 -5848±14 
11701±29 


3274±72 -6014±133 
11974±257 


Cij 


M c 1.00 0.971 
77 - 1.00 


1.00 0.962 
1.00 


1.00 0.963 
1.00 


1.00 0.967 
1.00 


1.00 0.960 
1.00 


Oi 


0.00776 0.00413 


0.00673 0.00363 


0.00735 0.00339 


0.00704 0.00364 


0.00628 0.00328 



TABLE V: Non-spinning binary fit, using an alternate technique: As Table ITTA except the effective fitting function 
allows Pmax to be a parameter. In this example, this more generic fit produces effective Fisher matrices that more closely 
correspond to conventional stationary-phase results for leading-order emission. 
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PN ° rd ° r 








le 


















h 


igher-orde 








iterative 




P ~ 1 






parameter 




M c V 


X 


M c 


V 


X 




M c 


V 


X 


M c 


n 


X 


M c 


V 


X 




M c 


3686 -1652 


-1007 


3567±22 


-1570±14 


-975. 6±5 


9 


2837 


-2423 


-669.6 


4217 


-2117 


-1177 


4129±34 


-2083±25 


-1157±9.3 




V 


1237 


515.5 




1170±9 


492. 4±3 


5 




2357 


612.3 




1806 


685.6 




1765±18 


670.5±6.4 




X 




283.8 






275.5 ±1 


6 






163.8 






340.0 






335.5±2.6 




M c 


1.00 -0.947 


0.996 


1.00 


-0.957 


0.997 




1.00 


-0.981 


0.994 c 


1.00 


-0.929 


0.994 


1.00 


-0.935 


0.995 




V 


1.00 


-0.969 




1.00 


-0.974 






1.00 


-0.995 




1.00 


-0.958 




1.00 


-0.962 




X 




1.00 






1.00 








1.00 






1.00 






1.00 






0.0291 0.0180 


0.135 


0.0323 


0.0201 


0.149 




0.0512 


0.0621 


0.440 


0.0227 


0.0131 


0.104 


0.0238 


0.0137 


0.108 



TABLE VI: Effective fitting parameters for an aligned-spin binary. As Table Hill for a BH-NS binary with x = 1- This 
binary has one new parameter (x). For all parameters shown here, a local quadratic form is a good approximation to the 
ambiguity function; this table provides the fitting parameters, treating the line of sight as fixed. For comparison purposes, this 
table provides both our standard simultaneous least-squares technique (applied to M c ,n, \, with all other parameters fixed) and 
an iterative technique applied to successive low-order subspaces. Table [VTT1 provides the results for a more generic fit including 
line-of-sight dependence. 



fitting scale 


P > 0.99 


parameter 




M c 


V 


X 


L 






M c 4388±54 


-2256±40 


-1230±15 


-0.2527±0.47 


-0.5585 ±0.38 


(Fy)eff 


'/ 




1869±28 


719.0±9.8 


0.3000±0.31 


1.122±0.23 


X 
i 

<$> 






356.2±4.2 


0.1098 ±0.14 
0.03279±0.0009 


0.3396±0.11 
-0.001972±0.0007 
0.02530±0.00025 




M c 


1.00 


-0.960 


0.997 


-0.206 


-0.622 




V 




1.00 


-0.976 


0.189 


0.567 


Cij 


X 

L 
<$> 






1.00 


-0.204 
1.00 


-0.616 
0.186 
1.00 


Oi 




0.0324 


0.0177 


0.147 


0.566 


0.818 



TABLE VII: Effective fitting parameters for an aligned-spin binary. Higher-order: As Table IIV[ but for a BH-NS 
binary with \ = 1. To a first approximation, the Fisher matrix roughly separates the intrinsic and extrinsic parameters. The 
weak covariances shown here are extremely susceptible to small changes in the Fisher matrix coefficients; only a single preferred 
value is shown. 



produce a dramatic qualitative change in how well the 
line of sight n can be measured. 

Both with and without higher harmonics [Fig. [5] , the 
line of sight is very difficult to measure, particularly 
at the expected relevant scale P ~ 0.99 (i.e., SNR of 
around 10). In both cases, the ambiguity function has 
a broad, extended, asymmetric extremum. In the ab- 
sence of higher harmonics, the ambiguity function cannot 
be usefully described by a locally quadratic approxima- 
tion, even a effective one. Nonetheless, by understand- 
ing the expected dependence on angle (and by adopt- 
ing coordinates in band), we can propose a physically- 
well-motivated fitting function, both for the purely an- 
gular dependence and for the correlations between line 
of sight and other parameters [Eq. (|30l) ]. This fitting 
function works extremely well when higher harmonics 
are neglected. When higher harmonics are included, a 
quadratic approximation sufficics, as the local extremum 
is much narrower. In the latter case, Table IIVI pro- 
vides the fitting parameters needed to reconstruct the 
full multidimensional fit. In fact, our well-chosen refer- 
ence frequency produces a nearly separable fit, with zero 



off-diagonal terms [e.g., r a jy — O,G0 C ~ in Eq. ([5TJ])]. 
We will return to this simple structure frequently below. 

Several effects besides higher harmonics also intro- 
duce fine-scale structure into the ambiguity function. As 
demonstrated in Fig. [5J the choice of reference frequency 
can introduce strong, scale-dependent features into the 
ambiguity function. We chose a reference frequency at 
100Hz to reduce its effect, but have not eliminated it com- 
pletely. The nonzero overlap between the (I, m) = (2, ±2) 
modes is another such effect. Hence, we are not surprised 
that our effective fitting parameters change as we reduce 
the range of P used in the fit, even in the absence of 
harmonics; see Table UITl 

The effect of scale dependence is fairly mild: the 
cigendirections for P > 0.99 and P > 0.999 agree, only 
the eigenvalue scale changes. For comparison, we also 
considered an alternate fitting technique that allowed the 
single best fit point to have P max ^ lj see Fig.[3J and Ta- 
ble El While this method leads to aesthetically pleasing 
results similar to analytic calculations, this fit system- 
atically underestimates P near the maximum-likelihood 
point and does not completely eliminate the trend to- 
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pN order 


leading-order 


parameter 
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M c 


3234 


-281.1 


-612.8 
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-1.359 


3.726 


-10.18 




'/ 
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929.9 


-3.089 


-171.7 
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-6.253 
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(L ij Jeff 
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132.0 
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0.3616 
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#nj 
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Qjl 
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-0.1021 




<\> 
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- 


- 
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0.3829 




M c 


1.00 


0.407 


0.730 


0.144 


-0.0729 
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- 




'/ 


- 


1.00 


0.742 


0.729 


-0.183 


0.109 


- 


Cij 


X 


- 


- 


1.00 


0.752 


-0.180 


-0.133 


- 






- 


- 


- 


1.00 


-0.191 


-0.106 


- 




0NJ 
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- 




1.00 


0.0126 


- 




QJL 
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1.00 


- 


o~i 




0.00592 


0.00573 


0.0433 


0.0212 


0.116 
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pN order 










;her-order 






parameter 
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139.3 
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-0.1706 
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( L tj Jeff 




- 


- 
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307.9 


0.02350 


1.998 


-7.544 






- 


- 


- 




0.8305 


0.003262 


0.006150 




Qjl 


- 


- 


- 






0.3510 


-0.1029 






- 


- 


- 






- 


0.3951 




M c 


1.00 


0.543 


0.778 


0.309 


-0.154 


-0.177 


- 




'/ 


- 


1.00 


0.815 


0.786 


-0.265 


-0.00409 
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X 






1.00 


0.815 


-0.270 


-0.230 






P 








1.00 


-0.271 


-0.204 
















1.00 


0.0381 






QJL 












1.00 




Gi 




0.00626 


0.00626 


0.0495 


0.0237 


0.115 


0.188 





TABLE VIII: Effective fit for precessing binary: Iterative method: The effective Fisher matrix needed to reproduce the 
calculated ambiguity function, derived from an iterative fit; see Table ITXl for a full 7-dimensional fit. Both of these matricies 
have comparable lists of eigenvalues (3625, 908, 54, 4.9, 0.74, 0.29, ~ 0) and (3844, 999, 57, 4.1, 0.75, 0.28, ~ 0) [bottom] and with 
the nearly zero eigenvalue roughly corresponding to the (f> direction. Since the eigenvectors span several orders of magnitude, 
the correlation coefficients of this poorly-conditioned matrix are extraordinarily sensitive to small changes in the coefficients and 
have not been provided. Unlike the non-precessing case, the geometric and intrinsic parameters do not completely separate, 
even for these well-adapted coordinates. The one-parameter uncertainties at shown are calculated by omitting the (nearly 
unmeasurable) (f> direction from the fit. 



wards different fitting parameters on the smallest scales. 
We henceforth adopt P max = 1. 



To facilitate approximate comparisons with prior work, 
Tables IIIII and IIVI provide one-dimensional standard de- 
viations a and correlation coefficients c^. Unless other- 
wise stated, these quantities are derived solely from the 
Fisher matrix fits from that same table. For example, in 
Table IIIII the "measurement errors" <Xj follow from in- 
verting the 2x2 matrix shown, while in Table IIVI they 
follow from inverting a full 4x4 matrix. 



B. Aligned spin 

Repeating our effective Fisher matrix calculation for 
an aligncd-spin BH-NS binary leads to results qualita- 
tively similar to the zero-spin case. As previously, higher 
harmonics provide little additional information about in- 
trinsic parameters, here M C ,T), X'i see Fig. [5] For exam- 
ple, looking at data along a fixed line of sight, Table IVII 
shows that, on a component-by-component basis (r) and 
overall (oi), the two ambiguity functions with and with- 
out higher harmonics resemble one another. More di- 
rectly, Fig. [S] shows three-dimensional contours of nearly 
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-420.2 ± 21 


-640.4±2.4 


916.4 ±3.6 
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0.00584 
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TABLE IX: Effective fit for precessing binary: Simultaneous method The effective Fisher matrix needed to reproduce 
the calculated ambiguity function, derived from a full 7-dimensional least squares; see Table fyTTTl for an alternative iterative fit. 
Both methods produce comparable results, with comparable eigenvalue distributions. The difference between the leading-order 
result and a model including higher harmonics is small, though usually significantly in excess of our fitting parameter error 
(e.g., several standard deviations) and of the systematic differences between the two fitting methods. As expected, higher 
harmonics lead to smaller parameter correlations 



constant P function of M c ,rj,x for both leading- 
order emission (blue) and higher harmonics (red) . While 
higher harmonics clearly do provide more information 
about intrinsic parameters - the red surface is nested 
inside the blue - the addition of higher harmonics only 
marginally improves our ability to measure the least-well- 
constrained combination of M c ,rj,x- Precisely as in the 
non-spinning case, however, higher harmonics provide 
more information about the line of sight. Despite our 
line of sight providing some sensitivity to a symmetry- 
breaking (2,1) harmonic, a waveform with harmonics en- 
codes roughly similar information as a waveform without 
harmonics. Tables IVII and IVHI provide the effective fit- 
ting parameters we used to reproduce their ambiguity 



function. As in the non-spinning case, we find an approx- 
imately separable fit, though less so than before [see Cij in 
Table IVII] . Even allowing for weak correlations between 
intrinsic parameters and the line of sight, the overall pa- 
rameter covariances o~i with harmonics [Table IVIIj arc 
nearly unchanged from a model with only leading-order 
emission [Table lVl] . 7 



For leading-order emission, measurements of orientation and in- 
trinsic parameters separate completely. The intrinsic parameter 
uncertainties computed in Table |VT1 are therefore identical to the 
uncertainties for the corresponding parameters derived from a 
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FIG. 5: Comparison of the ambiguity surfaces between 
the leading-order (top panel) and higher-order (bot- 
tom panel) waveforms for the non-spinning binaries. 

SX is defined by a difference from the fiducial value in the Ta- 
ble (TT] Higher-order waveforms change the ambiguity surface 
and break degeneracies related to the inclination and orbital 
phase so these parameters can be marginally observable at 
the scale of P > 0.99. 



The aligncd-spin results at \ = 1 cannot be easily com- 
pared with the corresponding zero-spin results (x = 0). 
On the one hand, the component-by-component Fisher 
matrix coefficients like T m will differ significantly, as the 
waveform phasing changes as a function of x an d hence 
so does ( x (d, 1 1 i>) 2 ') . On the other hand, the aligncd-spin 
results allow a new parameter (spin) that was treated as 
fixed for the x = case, with nontrivial coupling to the 
other intrinsic parameters. The one-parameter uncer- 
tainties a, are dramatically increased by including this 



full 5x5 Fisher matrix. 




FIG. 6: Ambiguity function for aligned-spin: Points in 
M c ,n,x with P £ [0.99,0.991], shown as blue (leading-order) 
and dark red (with higher harmonics). The two long, narrow 
surfaces (ambiguity ellipsoids) these points cover illustrates 
gravitational wave measurements can constrain one combina- 
tion of M c ,n,x tightly (e.g., M c ); one less so (e.g., n); and 
one almost not at all. The two surfaces nearly agree, with 
the most significant change being a slight reduction in the 
least-well-determined direction. The close agreement between 
these two surfaces shows higher harmonics provide fairly little 
additional information to break this degeneracy. 



previously-neglected systematic effect. 

Both at leading and higher-order, our effective fit to 
the ambiguity function is complicated by the wide range 
of scales in T, even for fixed line of sight. As is well-known 
from previous Fisher matrix calculations with aligncd- 
spins [40l |4lj . the ambiguity function in M c ,r],x has 
strong correlations, producing a narrow and extended 
extremum. For the specific example described by our ef- 
fective Fisher matrix, the M c ,r],x submatrix has eigen- 
values ~ 4200, 400, 30, describing a strong hierarchy of 
scales. For our purposes, Fig. [5] demonstrates our fidu- 
cial aligned-spin binary cannot be distinguished from bi- 
naries with spin x ^ 0-8: for each x in this range, suitable 
combinations of M c , r\ exist with high overlap. For these 
extremely extended ambiguity ellipsoids, a fit that repro- 
duces P(Ao,A) over the full range in x might require a 
more generic functional form than the one adopted so 
far: a quadratic with constant coefficients in the neigh- 
borhood of the fiducial binary. Effectively speaking, how- 
ever, these additional degrees of freedom add little infor- 
mation with considerable expense. We will explore more 
complicated effective dependence in a subsequent publi- 
cation. 

As in the zero-spin case, we find significant differences 
on the smallest scales in P, in a fashion that depends 
on the reference frequency. Given the number of dimen- 
sions, complex functional form, sensitivity to numerical 
implementation like the sampling rate, and less immedi- 
ate observational relevance, we defer a detailed discussion 
of fine-scale effects to a subsequent paper. 
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C. Precessing spin: Case 1 

For the first of our two fiducial processing binaries, we 
find higher harmonics provide little added information 
beyond the constraints produced in the non-precessing 
case. This unfortunate but expected result can be seen, 
for example, from the one-dimensional covariances in Ta- 
ble IIX[ from the effective Fisher matrix coefficients F^ ; 
or from their comparable sequences of eigenvalues. That 
said, even in the absence of higher harmonics, the ambi- 
guity function for a precessing binary has simpler struc- 
ture than the non-precessing result, with reduced cor- 
relations among the "intrinsic" parameters; a somewhat 
less extreme hierarchy of scales (i.e., eigenvalues) 8 ; and 
roughly speaking a more quadratic ambiguity function. 

In fact, for a precessing binary the previous clear sep- 
aration between "intrinsic" and "geometric" parameters 
breaks down. As each instant the opening angle /3 of 
the precession cone of L around J depends on the rela- 
tive magnitude of L and S = J — L, as well as on their 
(nearly conserved) misalignment angle L ■ S. The mag- 
nitudes of L and 5" are essentially intrinsic parameters, 
characterizing the binary masses and BH spin; therefore, 
we expect the precession cone opening angle (3 to be inti- 
mately correlated with the intrinsic parameters. At the 
same time, the precession cone opening angle must be in- 
timately connected to the "geometric" parameters that 
define the orientation of the binary at the reference fre- 
quency: 9^j,ajL,4>- Specifically, the orientation of the 
precession of L relative to the line of sight is character- 
ized by the two angles 9n.j (setting the orientation of J) 9 
and ajL (fixing the orientation of L along the precession 
cone at the reference frequency). The orbital phase <p 
at the reference frequency fixes the binary's geometry in 
band. This cross-coupling between intrinsic and extrinsic 
parameters has quantitative consequences for the Fisher 
matrix. Roughly speaking, the two new eigenvalues in- 
troduced into T by allowing precession (here with values 
~ 5,0.7, associated with the /3 and ajL parameters) can 
be expected to lie between the very large (3000, 900, 100) 
and very small (~ 0.3, ~ 0) eigenvalues associated with 
the manifestly intrinsic (M c ,r), x) and extrinsic (0nj,4>) 
parameters. 

For non-precessing binaries, the choice of a 100 Hz ref- 
erence frequency nearly separated intrinsic and extrin- 
sic parameters. A suitable choice of reference frequency 



8 We can always rescale our eigenvalues by rescaling our coordinate 
units. However, in these units all parameters have a prior range 
of order unity, and physical meaning. As a result, our eigenvalues 
also have meaning: the coordinate combinations corresponding 
to the smallest eigenvalues have minimal impact on the overlap 
and can be ignored. 

9 The other angle needed to specify the orientation of J is equiv- 
alent to a rotation around the line of sight, i.e. the polarization 
angle. As the complex overlap maximizes over this angle, it is 
explicitly removed as a parameter. 



may yet further reduce the off-diagonal terms in our effec- 
tive Fisher matrix. For the present coordinates, however, 
we cannot cleanly decompose parameters into "intrinsic" 
and "geometric" parameters. Table HXl shows correlation 
coefficients Cy calculated by omitting the (nearly immea- 
surable) <j) coordiantc in T; no obvious block-diagonal 
form occurs. 



D. Precessing spin: Case 2 

By contrast to the relatively simple ambiguity func- 
tions seen so far, our second set of binary parameters pro- 
duces a significantly more complicated ambiguity func- 
tion, particularly in geometric parameters. For example, 
Fig.[7]shows the ambiguity function versus Qnj, ctjL and 
((3jL,ctjL) all other parameters fixed, for "case 2". As in 
case 1, we have a highly symmetric binary starting with 
L, J, and the line of sight in the same plane at our ref- 
erence frequency. However, in this case we start with L 
parallel to the line of sight, rather than perpendicular to 
it. The ambiguity function shows extreme sensitivity to 
the initial conditions and highly nonquadratic behavior. 
These differences occur despite the considerable similar- 
ity between case 1 and case 2: the two are, to an excellent 
approximation, the same configurations, just slightly off- 
set in time. By contrast, the change in ambiguity versus 
M c , ?y, x is well-described by a quadratic form. 

This extreme scenario demonstrates that even an effec- 
tive Fisher matrix has limits: sometimes, a more generic 
functional form including higher-order correlations must 
be used on relevant scales. That said, highly non- 
quadratic behavior only occurred for a high-symmetry 
binary. We expect typical binary initial conditions will 
produce nearly quadratic ambiguity functions. 



V. CONCLUSIONS 

In this paper we have used case studies of a coher- 
ent, two-detector ambiguity function P(Xq, A) to estimate 
how much and what kind of information higher harmon- 
ics provide about BH-NS binaries. Given the high dimen- 
sion of and severe degeneracies that plague the problem, 
we perform a tractable, idealized calculation instead of 
the straightforward but less easily understood explicit 
source in a real multi-detector network. Specifically, we 
place a single source directly overhead an idealized detec- 
tor pair, equally sensitive to both polarizations. For all 
binaries, we find that higher harmonics provide little ad- 
ditional information about the binary's intrinsic param- 
eters. Instead, at best they provide information about 
the orientation of the source relative to the line of sight 
(e.g., i, (j) for non-precessing binaries). Notably, higher 
harmonics make it easier to exclude a non-precessing bi- 
nary with L ~ ±n. 

When possible, we estimate the two-point function 
P(Aq, A) on scales of interest (1 — P < 1/p 2 ) by a locally 
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FIG. 7: Symmetry can produce highly nonquadratic 
behavior: Plot of the ambiguity function when changing only 
On j and ajL (top panel) and /3jl,cjl (bottom panel), for 
the high-symmetry binary parameter set "case 2" . 



quadratic function 1 — V a p8\ a 8X^ /2, whose coefficients 
we denote the "effective Fisher matrix." Using coordi- 
nates A adapted to the most sensitive frequencies of our 
network, we show that to an excellent approximation the 
Fisher matrix for non-precessing binaries separates into 
intrinsic (masses and spins) and extrinsic (geometric) pa- 
rameters. For precessing binaries, the angles describing 
the precession cone are intimately connected with the in- 
trinsic masses and spins as well as the geometric param- 
eters. Therefore, they correlate with both intrinsic and 
extrinsic parameters and entangle the two previously dis- 
tinct sets of parameters. In future work we will compare 
our results to detailed posterior parameter distributions, 
as computed by Markov-chain Monte Carlo simulations. 
Our results suggest suitable coordinates will substantially 
simplify the interpretation of these posteriors. 



In addition to our main results, this paper also provides 
technical suggestions of broad interest to the data analy- 
sis community. First, rather than performing the con- 
ventional (single-polarization, real) overlap maximized 
over time and orbital phase, we use a coherent two- 
polarization overlap maximized over time and polariza- 
tion [42| . This overlap is far more discriminating than 
the single-polarization result and, lacking a maximum 
over the polar emission direction </>, is well adapted to 
asymmetric situations such as precession or the presence 
of higher harmonics. Second, we illustrate the impor- 
tance of choosing coordinates adapted to the network's 
sensitive band and to the binary's in-band geometry. We 
have provided concrete examples of the consequences of 
poor choices for coordinate conventions, which can lead 
to pathological behavior in the complex overlap and am- 
biguity function. Finally, we recommend adopting an 
effective Fisher matrix, derived by fitting the ambiguity 
function on the scale at which variations in the param- 
eters could be plausibly detected. We have shown that 
the conventional approach to the Fisher matrix can be 
sensitive to unobscrvablc fine-scale structure which can 
give misleading results. By contrast, simply by comput- 
ing the ambiguity function over the area of interest, then 
fitting, one can explicitly verify whether a quadratic ap- 
proximation even applies, as well as assess its error. 
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Appendix A: Comparison between the Real and 
Complex overlaps 

One detector can not be sensitive to both polariza- 
tions, as it measures a single (real) strain variable h. The 
conventional definition of the overlap [Eq. therefore 
provides an inner product for a single real data sequence. 

The complex overlap [Eq. @], however, uses infor- 
mation about both polarizations [33| . To illustrate the 
differences between these two diagnostics, Figs. |S] and 
|9] show the overlaps (ho\h) between our fiducial non- 
spinning waveform and an identical binary, except for 
one parameter. In both cases, these overlaps are maxi- 
mized over event time and polarization (for the complex 
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overlap) or orbital phase (for the real overlap). 

First and foremost, Fig. [5] shows that the real overlap 
does not change as inclination (l) or the orbital phase 
(</>) are varied. In other words, with one detector, we 
cannot identify the inclination, which measures the rel- 
ative amplitude of h + to h x . Also, by maximizing over 
<t>, we lose information about it. The conventional single- 
detector overlap therefore provides no information about 
how well we can measure these parameters in a network 
sensitive to two polarizations. 

By contrast, because the complex overlap explicitly 
uses two polarizations, it can identify the inclination; see 
the dotted line in Fig. [8l For example, because it can dis- 
tinguish between left- and right-handed sources, the over- 
lap between antipodal directions (i.e., t = 0, L = tt/2) is 
zero. Moreover, if the source orientation is different from 
0, then the waveform carries information about the or- 
bital phase, which the complex overlap easily identifies. 
As a concrete example, if i = tv/2 for the fiducial signal, 
the orbital phase (j> is measurable at high SNR. When 
l = 7r/4, <j) is unmcasurable for the leading-order wave- 
forms (for the analytic formula, see Eq. (Bl) in (HI) but 
marginally measurable for the higher-order waveforms 
(see Fig.©. 

On the other hand, when intrinsic parameters are 
varied, the complex and real overlaps largely agree, even 
though the complex overlap has more information avail- 
able; see Fig. HO 

For the figures and discussion in this paper, we com- 
pare the real and complex normalized overlaps, dividing 
each network response by the network signal to noise. 
For the two identical detectors used in the complex over- 
lap, the network SNR is just a quadrature sum of two 
detecters; see Eq. ([3]), 

p 2 = pI + pI. (Al) 

By contrast, for the real overlap, the relevant SNR is just 
p 2 + . For an identical source, a two-detector network has 
higher overall SNR. 

One point of this paper is to compare the real and com- 
plex overlaps from a parameter estimation point of view. 
Starting from Eq. (jlip . the real overlap enters directly 
into the expression for the posterior. As a result, con- 
tours of the real ambiguity function (ho\h) should closely 
correspond to contours of the posterior parameter distri- 
bution, for measurements limited to a single polarization 
and known sky location. By contrast, for a source with 
known sky location seen by a network with comparable 
sensitivity to both polarizations, the complex overlap en- 
ters directly into the expression for the posterior. A real 
network will have unequal sensitivity to two polariza- 
tions. We therefore expect the real posterior will resem- 
ble some average between the posteriors estimated using 
the single-detector and network overlaps. In a subsequent 
publication we will compare our results with posteriors 
computed by MCMC, to quantify how well our simple 
estimates do at characterizing measurement accuracy. 
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FIG. 8: Comparison of the ambiguity functions be- 
tween a real and complex overlaps for the extrinsic 
parameters. Dotted line is calculated by changing l where 
the fiducial value is and other parameters are the same as 
in Table [T] Dashed line is calculated by changing <j) where 
l = 7r/2 and other parameters are the same as in Table [J The 
ambiguity surface for the real overlap is flat. 
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FIG. 9: Comparison of the ambiguity functions be- 
tween a real and complex overlaps for the intrinsic 
parameters. Parameter values are summarized in Table [U 
Thick lines are calculated by changing Af c , others by changing 
V- 



Appendix B: Reference frequency and fine-scale 
structure in the ambiguity function and systematic 
errors 

Even for a zero-spin binary with only leading-order 
harmonics, the choice of reference frequency at which the 
parameters are defined significantly influences the struc- 
ture of the ambiguity function P(Ao, A). As concrete ex- 
amples, Figs. HO and [TT1 show that specifying the orbital 
phase at the start or end of the waveform introduces ad- 
ditional structure on physical scales into the ambiguity 
function for the most well- determined and physical pa- 
rameters, the chirp mass M c and mass ratio r/. Even more 
troubling, Fig. [TU] shows that an ill-chosen reference fre- 
quency can introduce extremely fine-scale structure into 
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the ambiguity function, for orientations away from z. By 
contrast, a reference frequency / re f close to the half-power 
point of the detector reduces these effects, where / re f is 
estimated by 



df 

/ref dp 2 



JM/)I S 



Si, 

D dl_ 



(Bl) 
(B2) 



for \h(f)\ oc /~ 7 / 6 the standard restricted amplitude. 

What introduces this severe dependence for such a 
vanilla waveform? The significant accumulation of or- 
bital and hence waveform phase between the coordinates' 
base point and the detector's sensitive band. Roughly 
speaking, the detector is sensitive to the configuration of 
the binary as it crosses through its sensitive band. This 
trajectory can be characterized by some instantaneous 
parameters A. By contrast, the waveform at a signifi- 
cantly earlier or later time has rotated roughly A$/7r 
times between that frequency and the observed one. As 
a result, derivatives of the waveform relative to A dif- 
fer from derivatives relative to A by a term of order 
hd5$/dA a . As a result, the more sensitive the waveform 
is to a parameter, the more severe the absolute impact of 
adopting a poor reference frequency. 

While we have chosen a single reference frequency, 
we have not optimized it, as we anticipate no single 



choice will work well for all parameters. The early part 
of the inspiral carries far more information about the 
chirp mass and mass ratio, as these both impact the 
total number of cycles. For example, to leading-order 
rjvfcM c — ((^m c ^) 2 ) where the average (} corresponds to 
averaging in signal power (i.e., averaging over a distribu- 
tion oc dp 2 /df). Evaluating the derivative and writing as 
an integral over frequency, we see the corresponding in- 
tegrand to dp 2 /df in the definition of Tm c m c is weighted 
by a polynomial in / that is heavily biased towards low 
frequencies: 



M C M C OC 



If 



-10/3 



lower powers] dp 2 /df (B3) 



By contrast, spin effects enter at v 3 oc / past leading- 
order in phase, suggesting the spin-spin Fisher matrix 
components are less severely biased towards low fre- 
quency; 



Tss oc / / 



-10/3+2 



dp 2 /df 



(B4) 



Other components of the Fisher matrix arise from 
other averages over frequency and have other break-even 
points. In other words, for each measurable quantity, a 
generally distinct epoch of the signal provides the most 
information. No one breakeven point works for all. 



[1] Abbott et al. (The LIGO Scientific Collaboration). Nu- [15 
clear Instruments and Methods in Physics Research A, 

517, 154 (2004). [16 

http : //www. virgo . inf n. it] 

Abbott et al. (The LIGO Scientific Collaboration). Phys. [17 
Rev. D 77, 062002 (2008). [18 
Abbott et al. (The LIGO Scientific Collaboration). Phys. 
Rev. D 79, 122001 (2009). [19 
A. Buonanno, Y. Chen, and M. Vallisneri, Phys. Rev. D 
67, 104025 (2003). [20 
A. Buonanno, Y. Chen, Y. Pan, and M. Vallisneri, Phys. 
Rev. D 70, 104003 (2004). [21 
T. Damour, A. Gopakumar, and B. R. Iyer, Phys. Rev. 
D 70, 064028 (2004). [22 
Y. Pan, A. Buonanno, Y. Chen, and M. Vallisneri, Phys. [23 
Rev. D 69, 104017 (2004). 

C. Konigsdorffer and A. Gopakumar, Phys. Rev. D 71, [24 
024039 (2005). 

A. Buonanno, Y. Chen, Y. Pan, H. Tagoshi, and M. Val- [25 
lisneri, Phys. Rev. D 72, 084027 (2005). 

C. Konigsdorffer and A. Gopakumar, Phys. Rev. D 73, [26 
124012 (2006). [27 
M. Tessmer and A. Gopakumar, MNRAS, 374, 721 
(2007). [28 
M. Hannam, S. Husa, B. Briigmann, and A. Gopakumar, 
Phys. Rev. D 78 (2008) [29 
I. Hinder, F. Herrmann, P. Laguna, and D. Shoemaker, [30 
arXiv:0806.1037. 



K. G. Arun, A. Buonanno, G. Faye, and E. Ochsner, 
Phys. Rev. D 79, 104023 (2009). 

A. Buonanno, B. R. Iyer, E. Ochsner, Y. Pan, and B. S. 

Sathyaprakash, Phys. Rev. D 80, 084043 (2009). 

F. D. Ryan, Phys. Rev. D 52, R3195 (1995). 

P. Jaranowski and A. Krolak, Living Rev. Relativity 15, 

4 (2012). 

E. K. Porter and N. J. Cornish, Phys. Rev. D 78, 064005 
(2008). 

K. G. Arun, B. R. Iyer, B. S. Sathyaprakash, and S. 
Sinha, Phys. Rev. D 75, 124002 (2007). 
C. Van den Broeck and A. S. Sengupta, Class. Quantum 
Grav. 24, 1089 (2007). 

J. S. Key and N. J. Kornish, arXiv:1006.3759 

R. N. Lang and S. A. Hughes, Phys. Rev. D 74, 122001 

(2006). 

A. Klein, P. Jetzer, and M. Sereno, Phys. Rev. D 80, 
064027 (2009). 

R. N. Lang, S. A. Hughes, and N. J. Cornish, Phys. Rev. 
D 84, 022002 (2011). 

M. Vallisneri, Phys. Rev. D 77, 042001 (2008). 

S. Vitale and M. Zanolin, Phys. Rev. D 82, 124065 

(2010) . 

M. Vallisneri, Physical Review Letters 107, 191104 

(2011) . 

L. Kidder, Phys. Rev. D 52, 821-847 (1995). 

C. M. Will and A. G. Wiseman, Phys. Rev. D 54, 4813 

(1996). 



19 




-0.003 - 0.002 











— loo Hz L t : 









0.000 

6M C 



FIG. 10: Examples for a dependence of the ambiguity function on the reference frequency. Left: leading-order 
waveforms, i — 7r/4, Middle: higher-order waveforms, i = 7r/4, Right: leading-order waveforms, l = 7r/2. 
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